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Abstract 

O 
O 

■ The cross sections for J /if) absorption by it and p mesons are studied in 

a meson-exchange model that includes not only pseudoscalar-pseudoscalar- 
vector-meson couplings but also three- vector-meson and four-point couplings. 
\ We find that they are much larger than in a previous study where only 

pseudoscalar-pseudoscalar-vector-meson couplings were considered. Includ- 
ing form factors at interaction vertices, the J /if) absorption cross sections a n ^ 
and apj/j are found to have values on the order of 7 mb and 3 mb, respectively. 
Their thermal averages in hadronic matter at temperature T = 150 MeV are, 
t^j- \ respectively, about 1 mb and 2 mb. 
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I. INTRODUCTION 



A dense partonic system, often called the quark-gluon plasma (QGP), is expected to be 
formed in heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC), which will soon 
start to operate at the Brookhaven National Laboratory. Of all experimental observables 
that are sensitive to the presence of the QGP, charmonium is among the most promising 
ones. In particular, the dissociation of charmoniums in QGP due to color screening would 
lead to a reduction of their production in relativistic heavy ion collisions. The suppression 
of charmonium production in these collisions has thus been proposed as a possible signature 
for the formation of QGP [|TJ. Extensive experimental and theoretical efforts have been 
devoted to study this phenomenon |2r||. However, available experimental data on J /if) 
suppression in colliding systems ranging from pA to S+U are consistent with the scenario 
that charmoniums are absorbed by target and projectile nucleons with a cross section of 
about 7 mb ||. Only in recent data from the Pb+Pb collision at = 158 GeV/c in 
the NA50 experiment at CERN |4]] is there a large additional J /if) suppression in high Et 
events, which requires the introduction of other absorption mechanisms. While there are 
suggestions that this anomalous suppression may be due to the formation of QGP 
other more conventional mechanisms based on J /if) absorption by comoving hadrons have 
also been proposed as a possible explanation [§,0] . Since the latter depends on the values of 
J/i[) absorption cross sections by hadrons, which are not known empirically, it is important 
to have better knowledge of the interactions between charmonium states and hadrons in 
order to understand the nature of the observed anomalous charmonium suppression. 
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Knowledge of J /if) absorption cross sections by hadrons is also useful in estimating the 
contribution of J /if} production from charm mesons in the hadronic matter formed in rela- 
tivistic heavy ion collisions. Since the charm meson to J/if> ratio in proton-proton collisions 
increases with energy, it has been shown that J/if> production from hadronic matter may 
not be negligible in heavy ion collisions at the Large Hadronic Collider energies fTI| : I^|. To 



use J / if) suppression as a signature for the formation of QGP in these collisions thus requires 
the understanding of both J / if) absorption and production in hadronic matter. 

Various approaches have been used in evaluating the charmonium absorption cross sec- 
tions by hadrons. In one approach, the quark-exchange model has been used. An ear- 
lier study based on this model by Martins, Blaschke, and Quack [|1^] has shown that 



the J /if) absorption cross section a w1 p by pions has a peak value of about 7 mb at 
E'kin = ^/s — m n — m^p ~ 0.8 GeV, but a recent study by Wong, Swanson, and Barnes 
jl4l gives a peak value of only a n ^ ~ 1 mb at the same -Ekin region. In the perturbative 



QCD approach, Kharzeev and Satz [[15] have studied the dissociation of charmonium bound 
states by energetic gluons inside hadrons. They have predicted that the dissociation cross 
section increases monotonously with E^ in and has a value of only about 0.1 mb around 
-Ekin ~ 0.8 GeV. In the third approach, meson-exchange models based on hadronic effec- 
tive Lagrangians have been used. Using pseudoscalar-pseudoscalar-vector-meson couplings 



(PPV couplings), Matinyan and Miiller |16| have found <r n ^ ~ 0.3 mb at E^ in = 0.8 GeV. In 



a more recent study, Haglin [|17] has included also the three-vector-meson couplings (VW 
couplings) and four-point couplings (or contact terms), and obtained much larger values 
of J /if) absorption cross sections. Large discrepancies in the magnitude of a n ^ (as well as 
a p ^) thus exist among the predictions from these three approaches, and further theoretical 



studies are needed. In the present study, we use a meson-exchange model as in Ref. |17J but 
treat differently the VVV and four-point couplings in the effective Lagrangian and also take 
into account the effect of form factors at interaction vertices. 

Our paper is organized as follows. In Sec. [TJ, we introduce the effective hadronic La- 
grangian that we use to obtain the relevant interactions among J / if) and hadrons. The cross 
sections for J /if) absorption by 7r and p mesons are then evaluated. The amplitudes for 
the coherent sum of individual diagrams are checked to ensure that the hadronic current is 



conserved in the limit of zero vector meson masses. We then show in Sec. Ill the numerical 



results for the cross sections and their dependence on the form factors at interaction vertices. 



In Sec. [TV], we compare our results with other models, and give more discussions on form 
factors and the effect due to the finite p meson width. A summary is also given in this 
section. In Appendix A, we discuss the determination of coupling constants based on the 
vector meson dominance model. More detailed comparisons with the approach used in Ref. 



T7| are given in Appendix B. 



II. J /if) ABSORPTION IN HADRONIC MATTER 

A. Effective Lagrangian 

The free Lagrangian for pseudoscalar and vector mesons in the limit of SU(4) invariance 
can be written as 
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£ = Tr (d^d»P) - \Ti (FlF» v 



(1) 



where = d^V u — c^V^, and P and V denote, respectively, the properly normalized 4x4 
pseudoscalar and vector meson matrices in SU(4) 
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To obtain the couplings between pseudoscalar mesons and vector mesons, we introduce the 
minimal substitution 



d,P - v,p = d,P- l j- [V,P] , 



The effective Lagrangian is then given by 

£ = C + ^Tr (>P [P\ tf\ + [P, V M 



(3) 
(4) 



-Tr 



pt T/t 
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+ f Tr (w [yj, v?] + h) + ^ (t^> n [yj, 

The hermit icity of P and V reduces this to 

£ = Co + ^Tr (d^P [P, V„\) - ^Tr ([P, ^ 



(5) 



+ igTi (W [V„ V v \) + ^Tr ([V„ V t 



(6) 



Since the SU(4) symmetry is explicitly broken by hadron masses, terms involving hadron 
masses are added to Eq.(H) using the experimentally determined values. 



B. Effective Lagrangians relevant for J ftp absorption 

Expanding the Lagrangian in Eq. (||) explicitly in terms of the pseudoscalar meson and 
vector meson matrices shown in Eq.(^j), we obtain the following interaction Lagrangians that 
are relevant for the study of J/ip absorption by 7r and p mesons: 
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£nDD* = WttDD* D*^T ■ (DS^ ~ B^DlTj + H.C. , 

C^dd = W^dd V (DdpD - d^DD) , 
C^d, = ig^ D * D * [V (d,D w D* u - D*»d,Dl) + (d^Dl - iP u dpD* u ) D* 

Ar^DD* = 9-xipDD* V (plfD + DtD*^J ■ 7T , 

C-pDD = WpDD (Drd^D - d^DfD^ 
£pi>DD = gpipDD ip^DfD ■ pi , 



(7) 
(8) 

(9) 
(10) 



pD*D* — ig P D*D* 



(d^rDt - D* u rd^Dl) ■ + (D*"t ■ d^p v - d^r ■ p v ) 



+ D^if-iTd^Dl-f-d^Dl 



(11) 



In the above, f are the Pauli matrices, and tt and p denote the pion and rho meson isospin 
triplets, respectively, while D = (D°, D + ) and D* = (D*°, D* + ) denote the pseudoscalar 
and vector charm meson doublets, respectively. We note that exact SU(4) symmetry would 
give the following relations among the coupling constants in the Lagrangian: 



g-irDD* — g P DD — gpD*D* 



g^DD* — gpfD*o* 



4 ) g-ipDD 



QtpD*D* - , 



gpnpDD 



2VQ 



(12) 



C. J/ip absorption cross sections 

The above effective Lagrangians allow us to study the following processes for J/ip ab- 
sorption by 7r and p mesons: 



nip -> D*D, nip -> DD*, pip -> DD, pip -> D*D*. 



(13) 



The corresponding diagrams for these processes, except the process nip — > DD*, which has 
the same cross section as the process nip — > D*D, are shown in Fig. |TJ. 

The full amplitude for the first process nip — > D*D, without isospin factors and before 
summing and averaging over external spins, is given by 



Mi = Mf e 2u e 3X = Y, M u e 2,e 3 A, 



(14) 
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FIG. 1. Diagrams for J/tp absorption processes: 1) Trip — > D*D; 2) — ► -D-D; and 3) 
— > D*D*. Diagrams for the process 7r?/> — ► -D-D* are similar to (la)-(lc) but with each particle 
replaced by its antiparticle. 
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Similarly, the full amplitude for the second process pip — ► DD is given by 



(15) 



with 



MZ = 9pDD9ipDD (pi - 2p 3 ) M ( 1 2 ) (Pi - P3 + P4)", 

\t m D j 



M Z = -9 p dd9^dd (-Pi + ^PaY 



u — m 



-Pi -Ps+PiT, 



D , 



M% = 9p4>dd <T ■ 



For the third process pip — > D*D*, the full amplitude is given by 
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In the above, p.,- denotes the momentum of particle j. We choose the convention that particles 
1 and 2 represent initial-state mesons while particles 3 and 4 represent final-state mesons 
on the left and right sides of the diagrams shown in Fig. [II respectively. The indices p, u, A, 
and uj denote the polarization components of external particles while the indices a and j3 
denote those of the exchanged mesons. 

After averaging (summing) over initial (final) spins and including isospin factors, the 
cross sections for the three processes are given by 
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(22) 



(23) 



is the squared momentum of initial-state mesons in the center-of-momentum (cm.) frame. 

D. Current conservation 

The effective Lagrangian in Eq. is generated by minimal substitution, which is 
equivalent to treating vector mesons as gauge particles. To preserve the gauge invariance in 
the limit of zero vector meson masses thus leads to both VVV and four-point couplings in 
the Lagrangian. The gauge invariance also results in current conservation; i.e., in the limit 
of zero vector meson masses, degenerate pseudoscalar meson masses, and SU(4) invariant 
coupling constants, one has 



Mt- Xl Pjx, = 



(24) 
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FIG. 2. J/ij) absorption cross section (without form factors) as a function of the cm. energy 
of initial-state mesons. The solid curve represents the total contribution from both irip — * D*D 
and nip — > D*D processes. 



where the index Xj denotes the external vector meson j in process n shown in Fig. [I]. This 
then requires, e.g., A4l x p 3 \ = and M^^Viv = 0. In Appendix B, we shall explicitly check 
that the amplitudes given in Eq. fll5|), as an example, satisfy the requirement of current 
conservation. 



III. NUMERICAL RESULTS 



The coupling constant g n DD* can be determined from the D* decay width |fL9|j , and this 
gives g-nDD* = 4.4. Using the vector meson dominance (VMD) model, we can determine 
other three-point coupling constants. As shown in Appendix A, their values are 

9 p dd = 9pD*D* = 2.52 , q^dd = 9^d*d* = 7.64 . (25) 

For the four-point coupling constants, there is no empirical information, and we thus use 
the SU(4) relations to determine their values in terms of the three-point coupling constants, 
i.e., 

Qir^DD* = 9ttDD*9iPDD, 9pipDD = 2 9pDD9ipDD-, 9p^D*D* = 9 'pD* D*9ipD* D* ■ (26) 

To obtain analytical expressions for the cross sections so that they can be directly included 
into a computer code for numerical calculations, we have used the software package FORM 
p0|| to contract all Lorentz indices in Eq. (E2|). 



A. Without form factors 



Figure [| shows the cross section of J / ip absorption by ix and p mesons as a function of 
the cm. energy i/i of the two initial-state mesons. The cross section a^, shown by the 
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FIG. 3. Thermal average of J/ip absorption cross section (without form factors) as a function 
of temperature T. 



solid curve, includes contributions from both 7T0 — ► DD* and 7T0 — > D*D, which have same 
cross sections. It is seen that the three J/ip absorption cross sections have very different 
energy dependence near the threshold energy, max(mi + m 2 , +1714). While <j n ^ increases 
monotonously with cm. energy, the cross section for the process pip — > DD decreases rapidly 
with cm. energy, and that for the process pip — *> D*D* changes little with cm. energy after 
an initial rapid increase near the threshold. 

The thermal average of these cross sections in a hadronic matter at temperature T is 
given by 



{av) 



dz [z 2 - fa + a 2 ) 2 } [z 2 - (oi - a 2 ) 2 } K^z) a(s = z 2 T 2 ) 
4afi^ 2 («i)a^ 2 («2) 



(27) 



where «j = m^/T (i — 1 to 4), z = max(ai+a2, 03+04), K^s are modified Bessel functions, 
and v is the relative velocity of initial-state particles in their collinear frame, i.e., 



m\m 2 



E 1 E 2 



(28) 



As shown in Fig. ||, (cr^v) increases with increasing temperature, but (cr^v) varies 
only moderately with temperature. The contribution of the process pip — > DD to (cr p ^v) 
is seen to decrease slightly with temperature. These features can be understood from the 
energy dependence of the cross sections shown in Fig. |2| and the difference in their kinematic 
thresholds (i.e., m.3 + m.4 — mi — m 2 ), which are about 0.64, —0.14, and 0.15 GeV for the 
processes nip — > D*D(DD*), pip — > DD, and pip — > D*D*, respectively. The process 
■Kip — > D*D(DD*) has the highest threshold, while the process pip —>■ DD is exothermic 
and thus has no threshold. With a pion in the initial state, the process Trip — > D*D(DD*) 
requires very energetic pions to overcome the high energy threshold and thus has a small 
thermal average (cr^^v) at low temperature. At higher temperature not only are there more 



S 



energetic pions that are able to overcome the kinematic threshold but also the cross section 
for the process irip — > D*D(DD*) increases with the cm. energy as shown in Fig. § 
(cr^v) thus increases strongly with temperature. For the process pip — > DD, on the other 
hand, its contribution to the thermal average (cr p ^v) decreases with temperature because 
with increasing temperature there are fewer rho mesons at low energy, which contribute the 
largest cross section. The contribution of the process pip — > D*D* to (cr^v) changes slowly 
with temperature as a result of both the small threshold and the fact that the cross section 
only weakly depends on the cm. energy. 

Compared with the results of Matinyan and Miiller [IB|, we see that the inclusion of 
the VVV and four-point couplings increases by an order of magnitude. For the process 
pip — > DD, the decrease of its cross section after including four-point couplings is due to their 
destructive interference with the PPV coupling terms. The process pip — > D*D* is entirely 
due to VVV and four-point couplings and is seen to have a much larger cross section than 
that for the process pip — > DD. As a result, our effective Lagrangian including the WV 
and four-point couplings also significantly increases a p ^. 



B. With form factors 

To take into account the composite nature of hadrons, form factors need to be introduced 
at interaction vertices. Unfortunately, there is no empirical information on form factors 
involving charmoniums and charm mesons. We thus take the form factors as the usual 
monopole form at the three-point t channel and u channel vertices, i.e., 

h = J^, (29) 

where A is a cutoff parameter, and q 2 is the squared three momentum transfer in the cm. 
frame, given by (pi — Ps) 2 m . an d (pi — P4)c. m for t and u channel processes, respectively. 
We assume that the form factor at four-point vertices has the form 

a-G&Xsf?)- (3o) 

where Ai and A 2 are the two different cutoff parameters at the three-point vertices present 
in the process with the same initial and final particles, and q 2 is the average value of the 
squared three momentum transfers in t and u channels: 



n 2 _ y-il -r VP1 H4; J c . m . _ 2 , 2 /oi \ 

H — 2 -«,c.m. +P/,c.m. ■ [ 6L ) 

For simplicity, we use the same value for all cutoff parameters, i.e., 

AttDD* = ApDD = Ap/)*D* = AipDD = A</,£)*D* = A, (32) 

and choose A as either 1 or 2 GeV to study the uncertainties due to form factors. 

Figure |] shows the cross section as a function of the cm. energy without and with form 
factors. It is seen that form factors strongly suppress the cross sections and thus cause large 
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FIG. 4. J/ip absorption cross section as a function of the cm. energy of initial-state mesons 
with and without form factors. 



uncertainties in their values. However, the J/ip absorption cross sections remain appreciable 
after including form factors at interaction vertices. The values for and a p ^ are roughly 
7 mb and 3 mb, respectively, and are comparable to those used in phenomenological studies 
of J /if) absorption by comoving hadrons in relativistic heavy ion collisions P,|IU,|2T|. The 



thermal average of J /if) absorption cross sections with and without form factors is shown in 
Fig. |5|. At the temperature of 150 MeV, for example, (cr^v) and {(Jp^v) are about 1 mb 
and 2 mb, respectively. 



IV. DISCUSSIONS AND SUMMARY 



In our study, the effective Lagrangian shown in Eq. @ is obtained from applying the 
minimal substitution of Eq. (|||) and Eq. (f|) to the free Lagrangian. The resulting PPV, 
VVV, and four-point (PPVV and WW) interaction Lagrangians in Eq. (|6|) are exactly 
the same as those in the chiral Lagrangian approach |22] []. They are, however, different 
from those used by Haglin ||17|| . The differences are shown in detail in Appendix B. 

Values of the J/ if) absorption cross sections by hadrons obtained in our model are compa- 
rable to those from Martins, Blaschke, and Quack [113], Haglin [I7|, and Wong, Swanson, and 
Barnes [Of, but are much larger than those from Kharzeev and Satz [15] and Matinyan and 



Muller [16|. As shown in Fig. and Fig. BL our results without form factors are much larger 



1 Before checking the identity of our effective Lagrangian to the corresponding ones in the chiral 
Lagrangian approach of Ref. [22], one should take notice of the different normalizations for the 
coupling constant g and the meson matrices as well as the following typos in that paper. The 



first part of Eq. (A2) in Ref. ]g2) should be U = exp 
dpU -igA^U + igUA^, and the first part of Eq. (A5) should be C 
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FIG. 5. Thermal-averaged cross section of J/if) absorption as a function of temperature T 
with and without form factors. 



than those from Matinyan and Miiller [Tb|] because the latter only included pseudoscalar- 
pseudoscalar- vector-meson couplings. As to the energy dependence of J /if) absorption cross 
sections, our results shown in Fig. ||] for the case without form factors are similar to those 



from Matinyan and Miiller [JIS] and Haglin [T7], which are also based on effective hadronic 
Lagrangians. Including form factors weakens the energy dependence of the absorption cross 
sections as shown in Fig. [|. However, the decrease of the absorption cross sections with 
energy is still not as fast as in quark-exchange models [13jl^|. This difference could be due 



to the fact that meson interactions in our effective hadronic Lagrangian approach involve 
derivative couplings, leading thus to a strong momentum dependence in the matrix ele- 
ments, while the nonrelativistic potential used in the quark-exchange model does not have 
an explicit momentum dependence. Including the relativistic corrections to the quark-quark 
potential will be useful for further studying the energy dependence of the J /if) absorption 
cross sections in the quark-exchange model. 

Form factors involving charm mesons introduce significant uncertainties to our model 
based on hadronic effective Lagrangians, because there is little experimental information 
available. Four-point vertices appear in all processes in our study. If all vector mesons 
are massless, it is possible to determine the form factor at a four-point vertex once form 



factors at three-point vertices are chosen [23]. This is achieved through gauge invariance 



by requiring current conservation for the total amplitude that includes the form factors. 
Since the uncertainty of form factors involving charm mesons is already large for three-point 
vertices and the gauge invariance is only exact when all vector mesons are massless, we 
choose not to follow this more involved approach. Instead, we show the uncertainties due to 
form factors by using two different values for the cutoff parameters. 

We have used the centroid value for the p meson mass in this study. Since the p meson 
width in vacuum is large (151 MeV), the threshold behavior of pip processes may change 
significantly with the p meson mass. E.g., a rho meson with mass below 630 MeV changes 
the process pip — > DD from exothermic to endothermic, and the energy dependence of the 
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cross section near the threshold may thus change from fast decreasing (the dashed curve) 
shown in Fig. || to fast rising (similar to the dot-dashed curve). On the other hand, a 
rho meson with mass above 920 MeV changes the process pip — *> D*D* from endothermic 
to exothermic. We thus expect that the final value of the J /if) absorption cross section by 
rho mesons will be different once the p meson width is considered. However, the p meson 
spectral function is further modified in the hadronic matter produced in heavy ion collisions 
P^jj2"5| , so the effects of the rho meson width on J/ip absorption in hadronic matter are more 
involved. We therefore leave the effect of the p meson width on charmonium absorption to 
a future study. 

Finally, vector mesons are treated as gauge particles in our approach. Since the SU(4) 
symmetry is not exact, it is not clear to what extent they can be treated as gauge particles. 
An alternative approach based on both chiral symmetry and heavy quark effective theory 
may be useful in understanding the meson-exchange model we have used here. 

In summary, we have studied the J/ip absorption cross sections by it and p mesons in 
a meson-exchange model that includes pseudoscalar-pseudoscalar-vector-meson couplings, 
three-vector-meson couplings, and four-point couplings. We find that these cross sections 
have much larger values than in a previous study, where only pseudoscalar-pseudoscalar- 
vector-meson couplings were considered. Including form factors at the interaction vertices, 
the values for and a pi p are on the order of 7 mb and 3 mb, respectively, and their thermal 
averages at the temperature of 150 MeV are roughly 1 mb and 2 mb, respectively. These 
values are comparable to those used in the phenomenological studies of J/ip absorption in 
relativistic heavy ion collisions. Our results thus suggest that the absorption of J/ip by 
comoving hadrons may play an important role in the observed suppression. 
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APPENDIX A 



In this appendix, we determine the values of the coupling constants within the framework 
of the VMD model. In the VMD model, the virtual photon in the process e~D + — > e~ D + 
is coupled to vector mesons p, u, and J/ip, which are then coupled to the charm meson. At 
zero momentum transfer, the following relation holds: 

£ 1V 9VD 2 +D ~ = e. (Al) 

In the above, jy is the photon-vector-meson mixing amplitude and can be determined 
from the vector meson partial decay width to e + e~, i.e., 

Tvee ~ 3<' (A2) 

with the fine structure constant a = e 2 /4ir. The relative signs of 7y's can be determined 
from the hadronic electromagnetic current expressed in terms of quark currents . Since 



the virtual photon sees the charge of charm quark in the charm meson through the ipDD 
coupling, we have the following relations: 

lip9i>D+D- _ 2 l P g P D+D- luguD+D- _ 1 / aq n 
2 ~~ o e ' 2 1 2 ~~ o e- v 6 

m^p 3 rrip 3 

Similarly, one has, from the process e~D° — > e~D°, 

2 ~~ o e > 2 1 2 ~~ ~o e - v A4 J 

m^p 3 rrip 3 

Using g p o+D- — ~9 p d°d° = 9 p dd, 9wD+d- — 9vD°d° = 9ujDD, and 9^d+d~ — 9^d°d° = 9i/jDD 
from isospin symmetry, we then have 

lip9^DD _ 2 J p g p DD , 1u9uDD _ 1 l P 9pDD , 1^9u,DD _ 2 

2 ~~ o e ' 9 ' 9 ~~ o e ' 9 1 9 ~~ ~o e - 

m^p 3 rrip 3 rrip 3 

From the above equations, we obtain the following coupling constants: 

9 p dd = 2.52 , g uDD = -2.84 , g^ DD = 7.64 . (A6) 
We note that in Ref. |Tj| the same VMD relations for 9 p dd and g^DD as our Eq. (|A5|) are 



used but slightly different values, i.e., g P DD = 2.8 and g^oo = 7.7, are obtained. 

Equations similar to Eq. ( |A5|) can be written for kaons and pions in order to obtain gvKK 
and gv-TTTT- The resulting coupling constants, multiplied by the corresponding prefactors in 
the following SU(4) relations, are given in the parentheses for comparison: 

<7 p7nr (5.04)= 2g pKK (5M) = 2g pDD (5.0A) = -y^ DD (9.36) . (A7) 

We note that \g p7Tn \ is 6.06 if it is determined from the p meson decay width to two pions. 
It is seen that the predicted values differ only slightly from the above SU(4) relation except 
the coupling constant g^DD- This may indicate a sizable uncertainty in the ipDD coupling. 
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APPENDIX B 



In this appendix, we discuss in detail the differences between our approach and that 
of Ref. [T7). In particular, we compare the general effective Lagrangian for interacting 



pseudoscalar and vector mesons and the specific interaction Lagrangians for J / ip scattering 
by pion and rho meson in the two approaches. We also examine the condition of current 
conservation for the amplitudes derived from the two approaches. 



A. General effective Lagrangian 



In both approaches, one starts from the same free Lagrangian of Eq. ([!]). But our 
matrices for P and V differ from those of Ref. jlTJ by a factor of l/y/2 as given recently 
in Ref. | 28fl . For the minimal substitutions given in Eq. (|3|) and Eq. (|j) for obtaining the 
interaction Lagrangians, the first one is the same in the two approaches but the second one 
is different. Instead of the factor g/2 in the last term of Eq. (|4]), Ref. [|17j uses g. As a 
result, the effective Lagrangian given by Eq. (2) in Ref. [IT] has the following form: 



Z int = igTz {PVdpP - d^PV^P) + -# 2 Tr (PV^P - PV^PV^) 
+ igTr (d»V v V v ] + [W, H 8,V U ) + s 2 Tr (VV [V„ V v \) 

V2igT±(& t V v [V ll ,V v ]) 



^Tr(^P[P,^])-^Tr([P,^] 



yTr([^,K] 



^Bl) 



This Lagrangian differs from ours in the three- (VVV) and four- (WW) vector meson 
couplings. Compared with our Eq. we find that our VVV and WW terms are a factor 
of 2 and 4 smaller than corresponding ones in Ref. |T7J , respectively. 

We note that the above differences in the effective Lagrangians used in the two approaches 
cannot be due to a redefinition of fields. To see this, we rescale the coupling g and the fields 
P and V by c g , cp, and cy, separately, then the relative ratios of the PW, PPW, VVV, 
and WW terms are given by 



2 222o 3 j a 2 ■ 

CgCpC V , c CpC v , 2c g c v , and 4c c 



(B2) 



It is obvious that these ratios cannot be changed to 1 simultaneously, as one needs c g cy = 1 
from the ratio of the first two terms but 2c q cy = 1 from the ratio of the last two terms. 



B. Lagrangians for the J ftp interaction with pions and rho mesons 



After expanding the general effective Lagrangian, we have Eq. (|7J) for the itDD* inter- 
action, which should be compared with the following one in Ref. JTT~ 



Ldd* = ^9' n DD' [DTiD^d^i - dpDTilfyri - H.c. 



(B3) 



where we have use g' to label the coupling constant in Ref. [|T7J. It is seen that our coupling 
constant g n DD* is a factor of 2 smaller. Apart from the possible difference due to the 
definition of D field, we have the same C^dd* as Ref. JTTj]. 
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Compared to our ipDD interaction Lagrangian in Eq. (j8|), Ref. has 
W = w'^dd V \DdpD - (dfj,D)D . 



(B4) 



Apart from a possible sign difference in the definition of g^r>D, both have the same C^dd- 
Instead of our ipD*D* interaction Lagrangian in Eq. (0), Ref. [ff7| has 



(B5) 



Besides a possible sign difference in the definition of g^o*D*, we have two more terms in 
J~-ipD*D*, which involve the derivative of the J/ip field, than in Ref. [17]. We note that the 
cyclic form of our C^d*d* yields the following factor for the three- vector meson vertex in a 
Feynman diagram: 



(Pi - V2) 1 g^u + (P2 - P3)p9in + (ps - Pi)»g~/n, 



(B6) 



which looks exactly like the structure of the three-gluon vertex in QCD. A similar difference 
appears between our approach and that of Ref. |I7J for the pD*D* interaction Lagrangian. 
We note that there may be typos in Ref. |L7| for 1^d*d* and I p d*d*, as D*D* and D*D* 
cannot be both scalar. 

For the nipDD* interaction Lagrangian, ours is given by Eq . (|10D while that in Ref. [17 

is 



which is non-Hermitian and thus likely contains a typo. 



(B7) 



C. Current conservation for nip — > D*D 

As pointed out in Sec. p ID[ in the limit of zero vector meson masses, degenerate pseu- 
doscalar meson masses, and SU(4) invariant coupling constants, all amplitudes for the three 
J/i/j absorption processes shown in Fig. [l| should satisfy the current conservation condition 
of Eq. (0). Here, we consider the process irip — > D*D as an example and explicitly check 
the current conservation condition in both our approach and that of Ref. |17j . For simplicity, 
we take all meson masses to be zero. 



1. Our approach 



The three amplitudes for the process irtp — > D*D are shown in Eq. (15). Multiplying 
them by p 2u and omitting the common factor — g^DD* gipDD> we obtain 



Mf aP2 



P2 ■ (Pi - P3 + Pi) 



-2 Pl +p 3 ) x = (2 Pl -p 3 ) 



(pi +PaT 



u 

(Pi + Pa) is 



gap 



P2 - Psfg uX + (-pi +P2+ PA) x g' 3v + (pi +P3- PA) v g 



P2v 



u 



-P3P2 +p!(-Pi+Pa) X -ug^ 



-s + t 
2u 



i A 



P3-P1- Pa 



Mf cV 2v = Pi- 



rn 
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Their sum is given by 



Mfp 2u = (-^) p x 3 =► 0. (B9) 

As indicated in the last step, it goes to zero when contracting with the external polarization 
e 3 \ of the charm vector meson. It is also straightforward to verify the current conservation 
condition for M" x p 3 \ as shown later in Eq. ( |B16| ). 



2. Approach of Ref. iZ 



Using the interaction Lagrangians of Ref. [I7| as explained in the above, we obtain the 



following corresponding amplitudes for Eqs. 



), and 



m la 



9nDD* 



a'^DD (~2P1 +P 3 ) 



m 



m 



uX 
16 



9-kDD 



(Pi ~P3 +P*Y, 

(pi-Pi)a(Pl-P4)p 



u — m 



D" 



9a(3 



ni 



X 



-p§9 vX + {-pi + PiYg pu + (pi + p 3 

~9nDD*9ipDD 9 ■ 



m 



,i/A 
lc 



(BIO) 



Taking g'^, D * D * = 9'ipDD as i n Eq. (5) of Ref. |L7j and omitting the common factor 
9ttDD*9^dd/^, we obtain 



m laP2u 

mwP2v 



i^Pi-Pz)\ 

{pi+paY 



u 

(Pi + P4)/9 



-9 a/3 



-Pz9 vX + ("Pi + Pi) x 9 liv + (Pi + Ps - Pa) v 9 



X n fiv 



P2v 



U 



P3P2 + P2 (-Pi + P4) A - ug 



px 



-s + t 
2u 



P3 - Pi - P4 



Their sum is thus 



m u 1 x p 2 v 



-s + t 
2u 



) P3 - 3p 2 



-3p x . 



(Bll) 



(B12) 



When contracting with the external polarization e 3 \ of the charm vector meson, the first 
term vanishes but the second remains as shown in the last step. Therefore, the current 



conservation condition is not satisfied in Ref. [17 



To understand the above results, we compare the amplitudes given in Eq. ( pil| ) against 
Eq. ( |B8|) and note the following two differences, (i) In mi& of Eq. QB11|) , the two terms 
~P2 9 uX + P2,9^ v involving the four- momentum of ip in the ipD*D* vertex are missing, (ii) 
mic in Eq. (|B11| ) is a factor of —2 larger than ours [also see the comment after Eq. (fB7Q ]. 
In calculating m vX p2v l the difference in (i) does not matter as it accidentally gives {—P2P2 + 



P2P2) = 0- The failure of satisfying the current conservation condition in Ref. [|17] is thus 
due to the difference in (ii). 
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45 



6 



5 - 



3.5 



Eq.(Bll) 

Eq.(Bl 1)— difference in (ii) 

Q O Eq.(Bll)-differences in (i)&(ii) = Eq.(B8) 



P 




4.0 



4.5 5.0 
Vs (GeV) 



FIG. 6. Comparison of J/ip absorption cross section by pions. The dashed curve represents 
the result of Eq. ( [Blip based on the Lagrangian from Ref. [[17]] . The solid curve is the result from 



Eq. (Bll) by dividing the amplitude m\ c by a factor of —2, which reproduces the result from Ref. 
[17]. The circled curve is the result from Eq. flB8|) based on our Lagrangian. 



To see this more clearly, we show in Fig. |6] by the dashed curve the cross section for J / ip 
absorption by pions obtained from the amplitudes given in Eq. ( |B11|) and using the coupling 
constants g'^ DD * = 8.8 and gL DD = g'^r>*D* = 7-7 given in Ref. fllTf . However, the results 
are different from that shown in Fig. 1 of Ref. fl7l, which is reproduced here by the solid 
curve in Fig. ||. We have found that to reproduce the results in Ref. [[Tj]] requires dividing 
the amplitude m\ c in Eq. ( [Bll|) by —2. 

We note that although the current conservation condition is satisfied for m\ x p2 V when 
mic in Eq. (|B11| ) is divided by a factor of —2, the current conservation condition for m^p^x 
remains violated. This is shown explicitly in the following. 

From the amplitudes shown in Eq. ( |B11| ), we have, after omitting the common factor 



m laP3X 



m lbP3X 



P3 



-2pi +p 3 ) 



t 

(pi +Pi) a 
u 

(Pi + Pi)p 
u 



(Pi - P3 + P±Y = (Pi - P3 + Pi) 1 



9aP 



-P^ X 



Pi+Pd x 9 pu + {pi+P3-p*Yg 



v J3\ 



P3 ■ (-pi +P4)g up + Ps(Pl 

P3- 



m \cP3\ 



m£p 3X /(-2 
Their sum is 



Pi) 



(Pi + Pi) 



P3A 
S - 



2u 



(Pi - Pi) 

(B13) 



mi A P3A 



(Pi + Pi) 



+ 



s-t 
2u 



(Pi - Pi) 



(B14) 



which does not go to zero when contracting with the external polarization ei v . 
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On the other hand, if we also add the two missing terms in 1^d*d* according to our Eq. 



i.e., adding —p29 vX J rP x 9 liv to the ipD*D* vertex in m^b in Eq.( [Bll[ ), we then have the 



following additional contribution to rn^pzx- 



(Pi + Pi) 



u 



(-P2P3 + P2 ■ pzd' 



(Pi + Pi) 



-2 



2u 



P3 



(B15) 



Combining the above two results gives 

'8-t 



2u 



(pi ~Pi- P3) 



s - t 
2u 



-Pi) 0. 



(B16) 



As shown in the last step, it vanishes when contracting with the external polarization e 2u - 
The results after eliminating both differences (i) and (ii), i.e., using our amplitudes given in 
Eq. (|B8| ), are shown by the circled curve in Fig. |6], which is only slightly different from our 
results shown in Fig. due to the different value for g^DD (7-7 vs 7.64). 
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